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Abstract 

We study the possibility of experimental verification of the spin=2 nature of the Kaluza- 
Klein (KK) graviton which is predicted to exist in the extra-dimensional Randal- Sundrum 
(RS) warped models. The couplings of these gravitons to the particles located on or near 
the TeV brane is the strongest as the overlap integral of their profiles in the extra-dimension 
is large. Among them are unphysical Higgses (W^ and Zl) and KK excitations of the 
Standard Model (SM) gauge bosons. We consider the possibility to confirm the spin-2 
nature of the first KK mode of the warped graviton (Gi) based on the angular distribution 
of the Z bozon in the graviton rest frame in the gg^ G\ — > W KK (Z KK )W(Z) — > WWZ, 
gg— ► G\ — ► ZZ and gg^ G\ — > Z KK Z — ► ZZH decay channels. Using Wigner D- 
matrix properties, we derive the relationship between the graviton spin, signal angular 
distribution peak value, and other theoretically calculable quantities. We then study the 
LHC signals for these decay modes and find that with 1000 fb _1 of data, spin of the RS 
graviton up to ~ 2 TeV may be confirmed in the pp -> W KK (Z KK )W(Z) -> WWZ -> 3 
leptons + jet + $x an d pp — > ZZ — ► 4 leptons decay modes. 
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I. INTRODUCTION 



With the upcoming start of the CERN LHC, our quest for the physics beyond SM 
is likely to yield positive results. On the theoretical side two of the most important 
questions to be answered are the Planck-weak hierarchy problem and the flavor 
puzzle of the SM. The Randall- Sundrum model with a warped extra dimension [l| 
is just about the only theoretical framework which simultaneously addresses both 
these questions making it a very compelling model of new physics. Perhaps the most 
distinctive feature of this scenario is the existence of KK gravitons with masses and 
couplings at the TeV scale which therefore should appear in experiment as widely 
separated resonances 

The original RS model as well as all of its extensions are based on a slice of 
Ad5*5 space. At the endpoints of this five- dimensional space (<f) = 0,tt), two branes 
are placed which are usually labeled as an ultraviolet (UV) Planck brane and an IR 
(TeV) brane; and the large hierarchy of scales is solved by a geometrical exponential 
factor. Postulating modest-sized 5 th dimension with radius R and curvature k the 
TeV/Planck ~ e -k-nR ra ^j G f sca i es can De numerically obtained by setting kR ~11. 
In the original RS model all SM fields were localized on the TeV brane. The only 
new particles in this model were KK gravitons with no SM gauge quantum numbers. 
Later, in addition to the KK gravitons, a bulk scalar field with a ^-dependent 
vacuum expectation value (VEV) was shown to generate a potential to stabilize 
the R modulus [jij]. However, this model leaves higher- dimensional operators in the 
5D effective field theory suppressed only by TeV scale which, in turn, generates 
unacceptably large contributions to flavor changing neutral current (FCNC) and 
observables related to the SM electroweak precision tests (EWPT). A natural way 
to avoid this problem, proposed by , is to allow SM fields to propagate 

in the extra dimension. In this scenario there are KK excitations of SM gauge and 
fermion fields in addition to those of the graviton. These states have masses in 
the TeV range and are localized near the TeV brane. The SM particles are the 
zero-modes of the 5D fields, and the profile of a SM fermion in the extra dimension 
depends on its 5D mass. By localizing light fermions near the Planck brane and 
heavier ones near the TeV brane, the contributions to the FCNC and EWPT are 
suppressed by scales 3> TeV. As a consequence, the KK graviton whose profile 
is peaked at the TeV brane will couple mostly to the top quark, Higgs (or, by 
equivalence theorem, to the longitudinal W and Z bosons) 0, [ly, U|, and KK 
excitations of the SM fields. 

Thus, the promising channels to observe RS gravitons are those where produced 
gravitons are decaying to fields localized near TeV brane. The search for the KK 
gravitons using its decays to the top quarks was performed in j^]. Signals from gravi- 
ton decay to Wl pair, which subsequently decay into pure leptonic or semileptonic 



final states, were considered in 12J. The 4-lepton signal through the decay to a pair 



of Zls was studied in [lOj. Reconstruction possibility of the Z's via their leptonic 



decays makes this a uniquely clean mode. 
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In this paper we would like to address the issue of confirmation of the spin-2 
nature of RS gravitons (for the most recent survey of methods of measuring the spin 
of new physics particles at the LHC, see 13J). The conventional way to measure 
the spin of a new particle involves reconstruction of its rest frame using its decay 
products and studying the angular distributions about the polarization axis. Along 
these lines, in ^ generic sample of a 100 tt events was produced for a spin-0, spin- 
1 and spin-2 resonances in an attempt to distinguish spin of the resonance based 
on the angular dependence of the cross-section. We, instead, concentrate on the 
gg^ Gi -> W KK (Z KK )W L (Z L ) -> ZZH,WWZ and gg-^ d -> Z L Z L channels 
where, in addition to the ZlZl channel considered in the literature before, we have 
one KK partner of the W or Z and one longitudinally polarized W or Z in the 
final state. Due to W KK and Z KK presence in these new channels, invariant mass 
of their decay products should show resonant behavior. Reconstruction of these 
intermediate KK gauge bosons will be important to reveal the internal structure of 
the RS model. We will see that Z boson signal angular distribution in the graviton 
rest frame for all these modes peaks at 90 degrees to collision axis. Performing 
angular analysis using Wigner D-matrix, we will derive the relationship between the 
graviton spin, angular distribution peak value, and other theoretically calculable 
quantities. As our method only requires to measure this peak value, where most of 
the signal events will be concentrated, we may optimistically achieve this goal with 
a relatively low sample of O(10) events. 



II. MODEL 



We closely follow the model discussed in [lOj and briefly review it here. As dis- 
cussed above, we allow SM fields to propagate in the extra-dimension and distribute 
fermions along it to generate observed mass spectrum without introducing addi- 
tional hierarchies in the fundamental 5D theory. SM particles are identified with 
zero-modes of 5D fields, and the profile of the fermion in the extra dimension de- 
pends on its 5D mass. As was shown before 

BBS 

, all fermion 5D masses are 0(1) 
parameters with the biggest one, among the SM quarks, being that of the top quark. 
To specify the model even further, the top quark is localized near the TeV brane 



and the right-handed isospin is gauged [14j . We consider tji being on the TeV brane 



(see discussion of the other possibilities in [10J], for example). At the end of the day, 
we are left with three parameters to be measured experimentally. We define them 
as c = k/Mpi (the ratio of the AdS curvature k to the Planck mass), fi = ke~ nkR 
monitors gauge KK masses with the first few being (2.45, 5.57, 8.7...) x/x, and finally 
the parameter v = m/k which defines where the lightest fermion with bulk mass m 
is localized. For the tn on the TeV brane, v tR ~ 0.5; and parameters c and \x will 
remain free in our analysis. 
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A. Low energy constraints on model parameters 



Before proceeding further, let us briefly review constraints placed on the warped 
extra-dimension model with custodial isospin symmetry 14| , which we adopt in this 
paper. In the resulting setup of 14|, KK mass scale as low as ~ 3 TeV is allowed by 
precision electroweak data. Regions of parameter space that successfully reproduce 
the fit to electroweak precision observables with KK excitations as light as ~ 3 
TeV were also studied in 15]. Phenomenological consequences of the observed BB- 
mixing were discussed in 16j . In the model of 14J tree level exchange of KK gluons 
gives the dominant contribution to the 55-mixing. In 16j , the CP- violating effects 
on the Bd system were shown to provide M^ on >3.7 TeV constraint at 68% CL. 



18 . In 17 



Phenomenological constraints from lepton-flavor- violations were discussed in 17, 
"anarchic" Randall-Sundrum model of flavor was studied, and the 
minimal allowed KK scale of ~ 3 TeV was found to be permitted for a few points 
in the natural RS parameter space; but models with custodial isospin can relax 



these constraints. In 18|, extensive analysis of B — >• K*l + l'~ modes was performed, 
concluding that only the B —> K*ee decay have sizable new physics effects. With SM 
contributions being suppressed, current experimental bounds were translated into 
the lepton bulk mass parameters. For the first KK gauge boson mass of 2-4 TeV, 
10-20% deviation from the SM results were found. Top quark flavor violations and 
B-factory signals were also studied in 19|, |20j, |21| . Finally, enhanced contributions 
to AS = 2 processes generated by beyond the SM operators with (V — A) <g) (V + 
A) structure, present in these frameworks, may impose additional constraints (22| . 
Without further flavor structure these contributions were expected to place a lower 



bound on the KK gluon mass of 0(8 TeV) |23|, |24|. However, most recent studies of 



the flavor constraints on the new physics mass scale find that the KK gluon mass 



should generically be heavier than about 21 TeV 25 



Relentless attempts to lower KK-mass scale further still flourish on the mar- 
ket. On this road, a number of other models were proposed trying to improve the 
pro spects to discover KK-particles at the LHC One of them is a model presented in 
[26j with a somewhat surprising claim that KK masses as low as 1 TeV are consis- 
tent with all current experimental constraints. An interesting variant of the warped 
extra dimension based on 5D minimal flavor violation was recently proposed in 27 



19, 28 with a 



The model allows to eliminate current RS flavor and CP problem 
KK scale as low as 2 TeV. Closing the list of examples, a volume-truncated version 
of the RS scenario called "Little Randall-Sundrum (LRS)" model was constructed in 
29l | . With the assumption of separate gauge and flavor dynamics, this setup allows 
to suppress a number of unwanted contributions to precision electroweak, Zbb, and 
flavor observables, compared with the corresponding RS case. 

Summarizing, we may say that KK gauge bosons with masses below 3 TeV 
(which would imply tuq > 4 TeV) would be difficult to have in current theoretical 
constructions. If this is the case, signals at the LHC, confirming the RS idea, would 



be extremely difficult to find; and studies conducted in lOj, |24| and later in this 
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paper support this unfavored future. However, in view of the above discussion, it 
also seems plausible that these models are still being developed; and, therefore, it 
is not inconceivable that explicit construction(s) will be found which will allow KK 
masses lower than 3 TeV without conflict with electroweak precision experiments 
and/or with flavor physics. This attitude was taken in |30| and we in this paper will 
also adopt this point of view. 



B. Couplings of KK gravitons 

After these brief remarks we can write the couplings relevant to our discussions 
here. Since the graviton couples to the energy- momentum tensor T^, coupling 
of the nth level KK graviton to the qth and mth level gauge bosons has the generic 
form: 

j- C qmn {qm) (n) 

Mpi K) 
where the magnitude of the C qmn coupling constants depends on the overlap of the 
particle wavefunctions in the extra-dimension. 

Analytic expressions for the coefficients C qmn with the flat zeroth mode gauge 
boson profile may be found in (2J and for the Wl and Zl on the TeV brane we need 
to replace them with delta functions. We present resulting couplings in Tabled along 
with partial decay widths for dominant decay channels for the lightest KK (n=l) 



graviton which will be the focus of our analysis; see also [lfj. The WlWl, Z L Z L 
and hh decay channels illustrate equivalence theorem once again, which is valid up 
to (mw,z I 'm-a) 2 where mc is the graviton mass. 

Let us briefly explain the result for the T{G n — »■ Wkk^l) decay mode as it 
involves off-diagonal elements of the energy-momentum tensor of the gauge fields. 



Tha gauge boson mass matrix is 17 



^ £ a mn A^A^\ (2) 

m,n=0 

and for the TeV brane Higgs scenario the off-diagonal elements ao m that describe 
the mixing of the zero and the mth KK mode a 0m = y2nkR. Thus, the off-diagonal 
elements of the energy-momentum tensor are given by: 

T^ (m ' 0) = ^ | l -F^FfJ - m^^kRA^Af } - (3) 
| F PM F (o) + F f) F irn) _ m 2 w V2^kR(A^A^ + Af A^)) . 

To the leading order in the my? /m^K, this results in the partial decay rate: 



T(G n -> W KK W L ) = ^< /( H^), 

480tt m 2 Gn m 2 KK m Gn 
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TABLE I: Couplings of the first level KK graviton to the SM fields. The tR is assumed to 
be localized on the TeV brane. Parameter mf is the mass of n=l graviton, xf = 3.83 is 
the first root of the first order Bessel function and e = e k7rR . N c = 3 is number of QCD 
colors. 



SM fields 




Partial decay widths for n=l graviton 


gg(gluons) 


€ 

2-KkR 


negligible 


w L w L 


e 


(cxf) 2 mf/4807r 


Z L Z L 


€ 


(cxf) 2 mf/960n 




€ 


N c {cxf ) 2 mf /320vr 


h h 


e 


(cxf) 2 mf/9607r 


W KK W L 


e 


390(cxf ) 2 mf /960vr • {m w /mf) 2 


Z KK Z L 


€ 


390(cxf ) 2 mf /960tt • (m z /mf ) 2 



where f(x) = r 2 + (6r 2 + 20r + Q)x 2 + 14(2 - r 2 )x 4 + (6r 2 - 20r + 6)x 6 + r 2 x 8 , 
r = \JlixkR « 8.4, and we neglected W boson mass in the phase-space consideration. 

In the class of models we are working with, mf ~ 1.5mf K for the mass of the 
lightest KK graviton and the gauge fields j3] which translates into /(2/3) ~ 173. 
As graviton mass changes from 1.5 to 3 TeV (which will be the typical range for the 
graviton mass we consider in this paper), for our numerical estimates we take the 
J Br(C7 1 (2.25TeV A ) -> W KK W L ) w J Br(G 1 (2.25TeV A ) -> w 1/2 x Br(G n -> 

The suppression in the coupling of the graviton to the gluons follows because 
the gauge boson has a flat wavefunction, and thus its couplings to the graviton is 
suppressed by the volume of the bulk trkR m 35. For the same reason, the decay of 
gravitons to transverse W and Z bosons, as well as photons, are suppressed by this 
volume factor. The masses of the KK gravitons are given by m n = x n \i where x n is 
n'th zero of the first order Bessel function. Notice that we do not need qqG coupling 
as it is Yukawa- suppressed, and graviton production is dominated by gluon fusion. 

In this model the total width of the graviton is found to be Tq = 
which is split between 6 dominant decay modes to WlWl, ZlZl, tntR, hh, ZkkZl, 
and WkkWl m the ratio 2:1:9:1:0.5:0.5. Taking c ~ 1, the total graviton width is 
~ 7% of its mass and is very close to the corresponding width for RS KK Z' in the 
same model 24 . 



III. GRAVITON SPIN MEASUREMENT 



Now we discuss the strategy to confirm the spin-2 nature of the first KK mode 
of the RS graviton in our channels. Out of five possible polarization states of the 
graviton gluons can produce only | JJ Z >= |2 ± 2 > and | JJ Z >= 1 20 > states due 
to two facts: gluons do not have longitudinal polarizations and the total angular 
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momentum has to be equal to J=2 (where we have chosen beam axis to be in 
the z-direction). Now, suppose that the two gauge bosons from graviton decay 
are produced at the polar angle 9. We rotate the gluons-produced graviton state 



specified by polarization tensor e^ v (JJz) by this angle [31 



e^(2J z ) = J2 D J% (0,^0)6^(24), (5) 



J'z 



where e^(JJ^) is the graviton state with the z-axis along the direction of the decay 
products, and D^ Jt (a, 9, 7) =< JJ' z \R{a, 9, i)\JJ z >= e'^d^j, (9)e- iJz ~* is the 
Wigner D-matrix. Independent Wigner small d-matrix elements for the spin-2 state 
are presented in Appendix [32|. Now we may easily derive the angular dependence 
of the helicity amplitudes for our channels. They follow from Eqj5]for the |2 ± 2 > 
graviton state which is produced by | H — > and | — h > gluons states: 

V (2 ± 2) = 4 2 2»^(20) + ^(5)^(21) + 41i(^)^(2 - 1). (6) 

Now just use Clebsch-Gordan decomposition of the e' (20) and e' (2 ± 1) 
states in terms of l<g)l final spin states to observe that, for example, helicity 
amplitude A[g(X 1 )g(X 2 ) -> Z(X 3 ))Z(X 4 )\ = A XlX2 x 3 x 4 for A + _ o ~ d^ (cos9) , 

(2) (2) 

A + _o_ ~ d 2 i (cos8), and A + _o+ ~ a^l^cosO), where we have used Z boson in the 
final state for concreteness. Notice that we have not included the ^±22(^)^1/(22) and 
d±2_ 2 (9)e'^ u (2 — 2) terms in Eqj6]as W and Z from graviton decay have longitudinal 
polarization and, thus, these terms cannot contribute. 

Why e^(20) graviton state does not contribute? This is again due to the fact 
that a gluon is massless. If you allow a gluon to have mass, you will obtain additional 
helicity amplitudes proportional to the mass of the gluon in agreement with the 
above angular analysis. For example, someone would find that A ++00 ~ d^(cos9), 
A ++0+ ~ d'^} 1 (cos9), etc. 

Inherent to our analysis is the assumption that a graviton is produced essentially 
at rest so that its decay products are mostly back to back. The requirement to find 
the graviton center of mass frame will limit possible decay channels for the gauge 
bosons as we will see later. If the rest frame cannot be reconstructed, we need to 
look for Lorentz invariant angular correlations which would encode information on 
the spin of the intermediate resonance. We do not pursue this approach here. 

Now we use the fact that the Wigner D-matrix elements D J mk (a,/3, 7) form a 
complete set of orthogonal functions of the Euler angles a, (3, 7 (we use symbols j 
and J for the total angular momentum quantum number interchangeably): 



2tt pn i>2tt q 2 



da I sinpdp d'j D 3 m , k ,(a, ^,7)* D J mk (a, ^,7) = . 5 m , m 5 k , k 5 fj , (7) 
Jo Jo a7 -h 1 



2j 



to determine the spin of the resonance state. 
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Taking into account SM background events, we observe that measured normal- 
ized angular distribution is related to the graviton spin in the following way: 

da _ E»=0,±1 Ci X [ d 2i( COS °)} 2 + °backgd{cOs6) 

adcosO 2j^T Si=o,±i 

@backgd 

where we used the normalization for the Wigner small d-matrix, Cj's are parton 
level cross section expansion coefficients convoluted over gluon PDF's, and we sum 
over three polarization states of the final state gauge bosons. As we will see for 
all the channels considered below, signal Z boson angular distribution peaks at 
pseudorapidity 77 = and, consequently, we will apply EqJHJat this point. 



IV. APPLICATIONS 



We estimated SM background with the aid of the COMPHEP package 33[ . For 
our graviton signal we used Mathematica program and partially cross-checked them 
with COMPHEP. CTEQ5M PDF's were used throughout (in their Mathematica 
distribution package 34| as well as intrinsically called by COMPHEP). 



A. ZZ decay channel 

To see the method at work, let us start with the simplest example of gg — ► 
G n — > Z L Z L — > e + e~, discussed in jlO] where the distinctive 4-lepton signal 



allows the reconstruction of all the masses of the particles in the decay chain. The 
dominant SM background for this purely leptonic mode is the pp — > ZZ+X, and the 
clean four-charged-lepton signal makes this mode a "golden" one. For this process 
we have: 



s 2 ({3 2 - 2)sin 2 6 [2 

2 = 

= A__ 00 = 0, (9) 



s ip — z sin u z 9.^,9 , jo) , 
A + _ 00 = A_ +m = W 2 ; = J -s 2 ((3 2 - 2)d%(cos6) 



where (3 = a/1 — 4M§ / s is the Z boson velocity. 

Neglecting the SM background for a moment and using EqJHl we find that: 



as d Jjdcos9 = 2 ^d^ (cose)d^>(cose) = \d%{cose)d%{cos9), (10) 

and, thus, the height of the peak in the normalized signal angular distribution as in 
Fig{Th is characteristic of the spin of the resonance. For our case, distribution peaks 
at 5/2 x (a/6/4) 2 =15/16 and is independent of the graviton mass. 



S 



Normalized signal angular distribution 



Normalized total angular distribution 




FIG. 1: (a) Normalized angular distribution for the signal cr(pp — ► G\ — > ZZ) cross- 
section due to the 1st KK graviton mode and (b) Normalized total (i.e. including SM 
background) angular distribution for the <r(pp — > ZZ) cross-section for mf = 1.5 TeV 
(solid) and mf = 3 TeV (dashed) integrated in the mf ± Fq ZZ invariant mass window 
with c= k/Mpi = 1. 



The irreducible SM background to theZZ final state is dominated by qq annihi- 
lation as gluon fusion proceeds via loop and, thus, interference of the KK graviton 
signal with SM background is negligible. The background cross-section exhibits for- 
ward/backward peaking due to t/u channel exchange while KK signal concentrates 
in the central rapidity region(lo|. 

On FigfjTj we show the total (signal plus background) cross-section integrated in 
the mf ir^ ZZ invariant mass window for two samples, 1.5 TeV and 3 TeV, graviton 
masses. For these mass values, we find SM background at cos6 l =0 as ~ 6% and 
~23% of the signal, respectively. We observe that for 1.5 TeV case peak value at 
cos#=0 changed to rs 0.5 due to the fact that for this graviton mass dbackgd ~ ^signal 
in mf ± To ZZ invariant mass window lOj and, thus peak value reduced by ~ 
half after normalization (see also Table [TTJ) . For the 3 TeV mass, the peak value 
~0.2 as Ubackgd ~ 3 x a signa i. We may also impose pseudorapidity t] cut to reduce 
background, keeping the signal (almost) unchanged. For example, for the t] < 2 cut 



considered in [10J, we find peak values as ~0.8 and ~0.5 for the 1.5 TeV and 3 TeV 
graviton masses respectively. Also, the background may be further reduced using 
lepton angular distribution to distinguish longitudinally polarized Z bosons from 
RS graviton decay from SM background 35 . 



Finally, using numerical results from Ref. lOj and including Z — > t + t~ channel 
not considered there, we obtain statistics presented in Table JTTl We assume 100% 
efficiency for our clean 4-lepton signal. Poisson statistics CL to observe at least one 
signal event will be appropriate description if the number of background events ^ 
10. We see that for 1.5 TeV and 3 TeV gravitons with 1000 fb _1 of data, we expect to 
have ~ 130 events and ~ 1 event respectively. This implies that higher luminosities 
are needed to reach 3 TeV graviton KK mass (for 3 ab _1 SLHC discussed in the 



community, see for example 36|). The reason for optimism on the issue of the 



detection of the r's from Z decay is that ~500 GeV energy r's have a decay length 
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TABLE II: Signal pp — ► ZZ — ► 4 leptons cross-section (in fb) for the = 1.5 TeV 
and rriQ = 3 TeV with the corresponding leading SM background. Numbers in brackets 
correspond to r\ < 2 cut case. For the low number of events, ^ 10, Poisson statistics is 
an appropriate description and the corresponding confidence level is, therefore, used. We 
assume 100% efficiency for our clean 4-lepton signal. 



1.5 TeV 


No cuts 


r\ < 2 cut 


# of events/1000 fb" 1 


S/B 


S/v^B 


Signal G -> ZZ -> 4 lept. 


0.13 


0.13 


130 


1.3(6.5) 


13(29) 


SM ZZ^ 4 lept. 


0.1 


0.02 


100(20) 






3 TeV 


No cuts 


T] < 2 cut 


# of events/1000 fir 1 


S/B 


CL 


Signal G -> ZZ -> 4 lept. 


0.001 


0.001 


1 


0.33(1.25) 


57% (54% ) 


SM ZZ^ 4 lept. 


0.003 


0.0008 


3(0.8) 







of I = 7TC m 20 mm and therefore might leave visible tracks in the detector 37 
B. W KK (Z KK )W L (Z L ) decay channels 



For our next examples we need to consider the matrix element for the gg— > 
G n — > W kk (Z kk )Wl{Zl) in the helicity basis. Working in the parton center of 
mass frame, the result is: 



M(g a g b -> W L W 



KK\ 



^— ' A l,2,3,4 



Al A2A3A4 



A, 



27ikRfi 2 s — (m^) 2 + iT^m^ 
where helicity amplitudes relevant for our process are: 



^+-00 



{V2irkR((r 2 — l) 2 m^ — r 4 s 2 ) — 4r 2 sm 2 v }sm 2 # 



2r 3 



'121 



A_ +0 _ = 
A — 00 = ^4++o 



yi^r {(V27rM- l)rr& + {V2^kR + l)r 2 sj (1 - cos6)sin6 
Jl}^ {{y/2irkR - l)m 2 w + (V^kR + l)r 2 s| (1 + cos6)sin6 



A. 



A. 



A. 







with r = mw/m^ K . To obtain M(g a g b — ► ZlZ kk ) replace mw with vriz ■ We have 
checked our results with Ref.[35| where the process gg^ G n — > ZZ was considered 
(which translates into r=l in our notation), and we confirmed them. 

After straightforward calculation we arrive at the parton level cross-section: 



dajgg -> W L W 
dcosO 



KK\ 



\M\ 2 

5127TS 



'1 - 



m 



KK\2 
W I - 



(13) 
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TABLE III: Branching ratios of W' and Z' in the TeV brane Higgs scenario. 



Decay modes 


w 


Z' 


W L H 


0.51 




w L w L 




0.35 


W L Z L 


0.40 




Z L H 




0.60 


tt 




0.05 


tb 


0.09 





where we have neglected the W boson mass in the phase-space consideration. 

As W and Z KK subsequently decay, we need to know their main decay chan- 
nels which we now turn our attention to. 

From now on, we generically call nth KK states of W and Z as W and Z' . We 
consider the simplified single bulk SU{2)l case and take (t, 6)z, to have close to a 
flat profile and t^ on the TeV brane as they together do the best in satisfying the 
combined flavor- changing neutral currents (FCNC) and precision constraints. After 
that, the decay widths for the leading channels of Z' and W are 24j : 



Y(W' -> tb) = g2sMmw ' T(Z' -> tt) = 9 sm( k v + K A) m z> 

T(W^W L H) = ^^L t nz'^Z L H) = g J^^ (14) 
T{W> W L Z L ) = 3lMcWm w ^ v{z , ^ WlWl) _ g 2 SM clK*m z , 



192vr v ' 1927T 



where in the TeV brane Higgs scenario k = g n / 9sm ~ V^nkR ~ 8.4 is the coupling 
strength of nth KK state relative to SM SU(2)l coupling, and c w (s w ) is the cosine 



sine) of the Weinberg mixing angle. Notice that Ref. 24| assumed Higgs as A 5 



38l . |39 | and, thus, the IR brane coupling enhancement is equal to V irkR there. 
Also, using the values of tiZ' overlap integrals for the fermion profiles specified 
above [24|, we obtain Ky w 1/4 — 5s^/3 and ka ~ — 1/4 — s^. The enhancement of 



SM coupling for the decay channels in the last two rows of Eq{14] follows from the 
fact that all the participating fields have a profile peaked near TeV brane compared 
to transverse zero-modes of W and Z, both having a flat profile in extra-dimension. 
Moreover, W — > tb decay channel is not enhanced since (t, 6)l fields have a close to 
a flat profile. Corresponding branching fractions implied by Eq{T4] are presented in 
Table ED 

We focus only on the 1st KK mode of W and Z as the effects of heavier KK 
modes are suppressed. On Figj2]we present the total resonant cross-section a(pp — > 
W L W KK ) ps a(pp — > Z L Z KK ) integrated in the mf ± r G W L W KK mass window. 
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Using branching ratios in Table II III total cross-sections after KK state decays may 
be easily obtained. 

Total o~(pp->W W KK ) [fb] cross-section 



1 r \ : 

. 1 r \. : 

s 0.01 , 

A \. 

a 0.001 r 

o' 

0.0001 r ; 

o.ooooi [ 

1500 2000 2500 3000 3500 4000 4500 5000 

M G (GeV) 

FIG. 2: (a) Total signal cross-section a(pp — > WlW kk ) ~ a(pp — > Z^Z KK ) as a function 
of the 1st KK graviton mode mass integrated in the mf ± r G WlW kk mass window and 
with c=l. 

Thus, the five possible final states are: WlWlZl (which may come from both 
Z L Z' and W L W intermediate states), tbW L , ttZ L , W L W L H, and Z L Z L H. In 
this paper we will concentrate on WlWlZl and Z^Z^H states. WlWlH final 
state faces the challenge to reconstruct efficiently the W mass from the W decay 
products. We will not consider tbPKi and ttZ^ final states as both of them may 
additionally be produced through an s-channel KK gluon exchange which couples 
strongly to the tt pair. 

ZZH decay channel 



We are now in a position to discuss the more complicated case of the pp — > G — > 
ZZ KK — > ZZH final state where we have three independent helicity amplitudes 
involved (see Eq ri2i) compared to the above ZZ case where only one independent 
helicity amplitude survived. We assume that both Z's decay leptonically so that 
both Z masses can be reconstructed. Then, as one of the Z bosons comes directly 
from the graviton decay, it will have a bigger energy than the other one. We again 
would like to know how angular distribution of this Z may help to determine the 
spin of resonance its emitted from. We consider the ideal situation of pure signal 
events first and then add background events (which will depend on the mass of the 
Higgs) later. 

Again writing da signa i/dcos6 = ^i^i x [<4i (cos6)] 2 and using C\ = C_i, we 
have: 
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do "signal 
(J signaldcosd 



Cp x [eg 



[cosO)} 2 + d x ([4 2 i(cos0)] 2 + [4? (co^)] 2 
^TT x (Co + 2d) 



(15) 



and, thus, the peak value occurs in this case at 



(3/8C0+C1/2) 



0.77. The normalized 



j(C +2Ci) 

angular distribution is shown on Figj3^ and again is independent of the mass of the 
graviton. Obviously, FigBh applies to WWH and WWZ CclSCS as well because the 
helicity amplitudes are the same. 



Normalized signal angular distribution 



Normalized total angular distribution 



CD 
CQ 

20.4 




3 <D 0.8 

W 

t o 




(a) 



(b) 



FIG. 3: (a) Normalized angular distribution for the signal cr(pp — > G\ — > ZZH) cross- 
section due to the 1st KK graviton mode and (b) Normalized total angular distribution 
for mf = 1.5 TeV (solid) and mf = 2 TeV (dashed) with cuts as in Eqll61 



Let us now consider the SM background for the ZZH case. It will depend on the 
leading decay mode(s) of the Higgs which, in turn, will depend on the mass of the 
Higgs boson. Important features can be highlighted by considering m# = 120 GeV 
case for which the leading Higgs decay mode is H — > bb. Due to the large Lorentz 
boost of the Higgs, we expect 2 b-jets to merge; and, thus, we conservatively require 
to have 4-leptons + 1 tagged b in the final state. 

We consider 1.5 TeV and 2 TeV sample graviton masses and impose the following 
cuts: 

m G = 1.5 TeV : \n ZtH \ < 2, m G -T G < M ZZH <m G + T G 
m G = 2TeV: \rj z , H \ < 2, m G - 2F G < M ZZH < m G + 2F G , (16) 

where for 2 TeV case we doubled the ZZH invariant mass window to increase the 
number of events. 

We use a b-tagging efficiency of 0.4 with a rejection factor for light jets (u, d, s, g) 
R=20 . We use a charm rejection factor R c = 5. In addition, we used 
BR(H -> bb)=0.7 and Ee )Al> r BR ( Z ~* ~ °- L A11 this results in the cross- 

sections presented in the second column of Table ITVl We find a clear signal above 
the background for 1.5 TeV case and 80% CL for 2 TeV case. Notice that we used 
the efficiency/rejection parameters optimized for low transverse momentum of the 
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TABLE IV: pp — > ZZH — ► 4 leptons 66 cross-section (in fb) for the signal with mg = 1.5 
TeV and = 2 TeV and the corresponding leading SM backgrounds with cuts as in Eqfl6l 
and efficiency /rejection factors as discussed in the text. For the low number of events, < 
10, Poisson statistics is an appropriate description and the corresponding confidence level 
is, therefore, used. 



1.5 TeV 


Chits and b-taer 


# of events /500 ftr 1 


S/B 


CL 


Signal G -> ZZH -» 4 leptons +66 


0.0196 


9.82 


6.7 


99.9% 


SM ZZb^ 4 leptons + b 


1.6 xlO" 4 


0.08 






SM ZZqi — ► 4 leptons+ % 


1.5xl0~ 3 


0.75 






SM ZZg^ 4 leptons + g 


1.2xl0~ 3 


0.6 






SM ZZc^ 4 leptons + c 


5.6xl0" 5 


0.028 






2 TeV 


Cuts and b-tag 


# of events/1000 fb" 1 


S/B 


CL 


Signal G -» ZZH -» 4 leptons +66 


1.82xl0 -3 


1.82 


1.36 


61% 


SM ZZb^ 4 leptons + b 


6.65xl0" 5 


6.65xl0- 2 






SM ZZ(/£ — > 4 leptons+ 


5.24xl0~ 4 


0.52 






SM ZZg^ 4 leptons + g 


7.28xl0~ 4 


0.73 






SM ZZc^ 4 leptons + c 


2.41xl0~ 5 


2.41xl0" 2 







b-quark Pxb, and rejection is expected to improve for high Pxb which is the case at 
hand. Also, efficient reconstruction of Z mass from hadronic Z decay will increase 
the number of signal events as those modes have a bigger BR. 

On Figj3)D. we show the normalized angular distributions for 1.5 TeV and 2 
TeV graviton masses considered. We observe that for 1.5 TeV mass 0.77 peak value 
remains (almost) unchanged as the distribution is dominated by signal events, while 
for the 2 TeV mass value peak is less distinct. 

WWZ decay channel 



As discussed above, signal angular distribution is the same as in Figj3^ for 
this case because the helicity amplitudes are the same. Additionally, WZ or WW 
invariant mass presumably should have resonant (due to W' or Z') distribution; but 
we don't impose cuts on WZ or WW mass as we would like to keep our analysis 
as general as possible. At this point we have to decide on the decay modes of W 
and Z boson. We again allow Z decay leptonically so that we reconstruct Z mass 
efficiently, and we use the angular distribution of this Z for determination of the 
spin of the graviton. Now, if we allow both W's decay hadronicaliy, due to the huge 
Lorentz boosts of these W's, we pick up 2 leptons + 2 jets as a background for our 
decay mode which we find to be overwhelmingly bigger than our signal. Thus, we 
use (W — > jet) (W— > leptons) and (Z— > leptons) as our final state. 

For the leptonic W decay, due to small angular separation between missing 
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neutrino and charged lepton, we may estimate longitudinal (L) component of the 
v's momentum as: 



- (17) 

Pt, 

Using this collinear approximation, the momentum of the leptonic W is recon- 
structed and, thus, we can calculate the (presumably) resonant invariant mass of 
the WW or WZ system. In doing so, we assumed that leptons are coming from the 
W decay as the reconstructed leptonic W mass will be zero in the collinear approx- 
imation. Also notice that in this approximation, the Mww measurement error for 
the TeV energy W bosons is ~ mw/Ew ~ 0.1. 

We again consider 1.5 TeV and 2 TeV sample graviton masses and impose the 
following cuts: 

\r) Z ,w\ < 1, m G -T G < M wwz < m G + T G . (18) 

We have to remember that EqJH] is valid only if we integrate over whole angular 
coverage of the detector. Fortunately, as our signal concentrates in central rapidity 
region, even such a hard pseudorapidity cut changed the signal cross-section for both 
graviton masses only by about 8% which is in the range of experimental uncertain- 
ties. This cut also changed the peak value in the normalized angular distribution of 
FigBh, from 0.77 to 0.83 value. 



In addition, we use the result of 24J which finds that jet mass cut: 



65 GeV < M jet < 115 GeV (19) 

achieves acceptance fraction of 0.78 for the signal and 0.3 for the background events. 
TablefVl shows our results after all this cuts are imposed. 

Finally, on Figj4]we show the normalized angular distributions for the 1.5 TeV 
and the 2 TeV graviton masses considered. We observe that for the 1.5 TeV mass 
peak value of 0.83 for the case of zero background changed to about 0.7. This 
value can also be obtained applying EqJH] and using the fact that for this mass 
Cbackgd ~ & signal as can be seen in TablefVl and the fact that cr signa i(cos9 = 0) ~ 
1.68abackgd(cos9 = 0). For the 2 TeV mass value, peak is no longer seen due to the 
dominance of the background. 



V. CONCLUSIONS 



In this work, we have extended earlier studies of the discovery potential of 
warped gravitons at the LHC which concentrated on the gravitons decaying into 
the "gold-plated" Z L Z L channel, WlWl, channel and into the it pair. We have 
considered resonant production of the first RS KK graviton mode via gluon-fusion 
process followed by its subsequent decay to W KK \Z kk )Wl{Z l ) and Z L Z L pairs. 
We focused on confirmation of the unique spin-2 nature of the graviton using Z boson 
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TABLE V: pp WWZ 3 leptons + jet + $ T cross-section (in fb) for the signal with 
vtlq = 1.5 TeV and tuq = 2 TeV and the corresponding leading SM backgrounds with cuts 
as in EqlTSjand Eq[T2]and efficiency /rejection factors as discussed in the text. 



1.5 TeV 


Cuts 


# of events/300 fb" 1 


S/B 


S/y/B 


Signal G WWZ -> 3 leptons + jet + $ T 


0.10 


30 


1.16 


5.9 


SM WWZ^ 3 leptons + jet + $ T 


0.0026 


0.78 






SM WZg -> 3 leptons + jet + $ T 


0.0656 


19.7 






SM WZg^ 3 leptons + jet + $ T 


0.018 


5.4 






2 TeV 


Cuts 


# of events/1000 fb" 1 


S/B 


s/Vb 


Signal G -► WWZ -> 3 leptons + jet + $ T 


0.008 


8 


0.26 


1.44 


SM WWZ^ 3 leptons + jet + $ T 


6.8xl0~ 4 


0.68 






SM WZq -> 3 leptons + jet + $ T 


0.023 


23 






SM WZg^ 3 leptons + jet + $ T 


0.0072 


7.2 







Normalized total angular distribution 




FIG. 4: Normalized total angular distribution for a(pp — > WWZ) cross-section for mf 
1.5 TeV (solid) and mf = 2 TeV (dashed) with cuts as in EqUSJand EqlH 



angular distribution in the graviton rest frame for all these modes. We performed 
angular analysis using Wigner D-matrix in order to derive the relationship between 
the graviton spin, angular distribution peak value, and other theoretically calculable 
quantities. As our method only requires to measure this peak value, where most of 
the signal events will be concentrated, it may be possible to achieve this goal with a 
relatively low sample of O(10) events. In any case, our main aim in this work is to 
illustrate how our method can work, at least in principle. Using statistical results 
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of lOj for pp — > ZZ — > 4 leptons and our analysis of pp —* W KK W — *> WWZ — > 
3 leptons + jet + decay modes, we showed that with 1000 fb" 1 of data, these 
channels allow this number of events to accumulate for the RS graviton up to ~ 2 
TeV. As a byproduct of our analysis, we found that W KK \Z^ kk ^)Wl{Zl) graviton 
decay modes, which have not been studied before, have a Br comparable to the 
zero mode decay channels as summarized in Tables JltiVl These decay modes open 
new channels to search for the RS gravitons. As an extra bonus, reconstruction 
of intermediate KK gauge bosons in these modes will be important to reveal the 
detailed workings of the RS model. 
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APPENDIX A: SPIN-2 WIGNER SMALL D-MATRIX 



4? 03) 



cos 



4 

sin z (3, 



cos (5 



stn(3, 4 2 i(/3) 



1 — cos/3 



sinf3 



1 + cos/3 



3cos 2 p 



(2cos0-l), 



(1 - cos(3f 



$w = - 

^-(2cos(3 + l), 



sin/3 cos(3 
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